Bilayer Graphene Interferometry: Phase Jump and Wave Collimation 
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We theoretically study the phase of the reflection amplitude of an electron (massive Dirac fermion) 
at a lateral potential step in Bernal-stacked bilayer graphene. The phase shows anomalous jump 
of 7r, as the electron incidence angle (relative to the normal direction to the step) varies to pass 
±7r/4. The jump is attributed to the Berry phase associated with the pseudospin-1/2 of the electron. 
This Berry-phase effect is robust against the band gap opening due to the external electric gates 
generating the step. We propose an interferometry setup in which collimated waves can be generated 
and tuned. By using the setup, one can identify both the n jump and the collimation angle. 

PACS numbers: 73.23.-b, 73.63.-b, 81.05. Tp 



Introduction. — Monolayer graphene has attracted 
much attention due to potential application of its unusual 
properties [H. Its low-energy quasiparticles are massless 
Dirac fermions. They show Klein tunneling anoma- 
lous quantum Hall effects 0, 0|, electron optics behav- 
ior such as focusing [j| and collimation @, 0, H, @], etc. 
The Klein tunneling is understood by the chiralityof the 
quasiparticles Q , or equivalently by Berry phase [HI EH • 
There have been experimental efforts to observe it in a 
bipolar junction 0,111 and in an interferometry [lH ]. 

Bilayer graphene has properties quite different from 
the monolayer. For example, in Bernal-stacked bilayer, 
low-energy quasiparticles are massive Dirac fermions [15| . 
There have been studies on bilayer-graphene Klein ef- 
fects _[2| jQuantum Hall effects [16(, band-gap engineer- 
ing [ijj O EM EiJ 01 j etc. However, more studies on 
the bilayer may be necessary for the understanding of 
the features of the massive Dirac fermions such as their 
transport related with Berry phase 16|, 21 1 and collima- 
tion [2j. 

In this Letter, we theoretically study the reflection of a 
low-energy electron at a lateral potential step of monopo- 
lar (p-p or n-n) type in Bernal-stacked bilayer graphene 
[Fig. [T][a)]. The phase of the reflection amplitude is found 
to show anomalous behavior, as the electron incidence 
angle Q\ (relative to the normal direction to the step) 
varies. It shows an abrupt jump of 7r at 9\ = ±7r/4 when 
the step height is much smaller than the kinetic energy of 
the electron. Based on a reversal symmetry [Eq. (|2J)], wc 
attribute the 7r jump to the Berry phase associated with 
the pseudospin-1/2 of the electron. The jump becomes 
gradual, as the step height increases, due to the evanes- 
cent waves existing at the step boundary and breaking 
the reversal symmetry. We show that the phase jump can 
be detected in an interferometry setup [Fig. Hfb)]. We 
remark that the phase jump is robust against the band 
gap opening due to the electric gates generating the step, 
and that the setup does not require an external magnetic 
field, contrary to previous works @, 0, U , 14 , HI for the 
detection of Berry phase effects in graphene. 

Potential step. — In Bernal-stacked bilayer graphene, a 
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FIG. 1: (Color online) (a) Lateral monopolar potential step 
in Bernal-stacked bilayer graphene (xy plane). The arrows 
represent the propagation of electron waves, the dashed lines 
show the boundaries between the regions with different po- 
tential strengths Vi's, 61(2) is the propagation angle of the 
incident (transmitted) wave, and n's (ti's) are the reflection 
(transmission) amplitudes. As 9i varies, the phase of n can 
exhibit an abrupt jump of tt at 6\ — ±7r/4, a Berry-phase ef- 
fect, (b) Interferometry setup for the detection of the n phase 
jump of n, based on the interference between the two paths 
drawn in x 6 [0, d s ]. To see the incidence-angle (here, 8 C ) 
dependence of n, collimated waves are generated from a col- 
lection of waves with different incidence angle unc , by using 
the resonant filtering due to the barrier V c . The collimation 
angle C can be also identified, using the interference. 



low-energy electron states y) with energy E in the 

K valley is governed by Hamiltonian [l5[ Hk^k = E^!k, 



-a K -q + V(x), q = (-p\ +p 2 y , -2p x p y ).{\) 



The two pseudospin components of vI/k describe the lat- 
tice sites Ai and B2 of the bilayer, where A; and B; de- 
note the two basis sites of layer I = 1,2. &k = (c x , a y ) is 
the pseudospin operator for the K valley, p — (p x ,p y ) is 
the momentum relative to the valley center, 7 ps 0.39 eV 
is the interlayer coupling, v ~ 10 6 m/s, and a x<y z are the 
Pauli matrices. H K is valid for 0.0027 < \E - V\ < 7, 
where the trigonal warping [TEj is ignored. We study the 
K valley only, as the K' valley shows the same result, and 
ignore the intervalley mixing as V(x) slowly varies on the 
scales of the lattice constant and hv/j. Band gap due 
to the gates generating V(x) will be considered later. 
It is worthwhile to see a symmetry of Hk, which has 
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FIG. 2: (Color online) Dependence of the reflection phase 
6 ri and probability |ri| 2 (scaled by a 2 ) on 6\ for a = —0.01 
(solid), -0.1 (dashed), and —0.5 (dotted). Around <?i = tt/4, 
ri shows a jump of 7T, accompanied by |ri| 2 = 0. 



not been discussed in literatures. Hk is invariant under 
the antiunitary operator 6 (for O.OO27 < \E — V\ <C 7), 



— ia y TZ^±^/ 2 C, 



(2) 



when p is real (which is achieved when V\ — V% — ► in 
Fig.Q]). Here, TZ^± 7r / 2 is the operator rotating p by angle 
±7r/2 [i.e., {p x ,Py) — > (±P y , TPk) and thus <f — > -g], i(x y C 
reverses pseudospin, and C is the complex conjugation 
operator. We call this invariance as reversal symmetry 
hereafter, in the sense that O is exactly the same as the 
time reversal operator defined for a single valley [22| if q is 
replaced by p in Eq. (JTJ) . The symmetry and the Berry 
phase 7T [231 ] . associated with a loop encircling once the 
origin in the q space, can give rise to an interesting effect 
in a monopolar potential step (see below). This effect 
corresponds to the Klein tunneling in bipolar monolayer 
graphene, which results from the time reversal symmetry 
and the Berry phase tt due to a p-space loop [HI EH ■ 

Now, we consider a potential step V(x) of monopolar 
type, where V(x) = V\ for x < and V2 for x > 0, 
and study the reflection amplitude r% of a plane wave 
incoming to the step with energy E and incidence angle 
#1 [Fig-Ufa)], where E > V u V 2 (E < Vi, V 2 ) for a n-n (p- 
p) step. Due to the translational invariance along y, the 



wave is described by >Fk. 



v i>K,k y (x). 4>K,k v (x) l is 



a superposition of propagating and evanescent waves [2| , 



yl + sin 2 6<i + sisin6>i 
Si\/ 1 + sin 2 #i — sin6*i 



for x < 0, 



tie 



-s 2 e 



+ 6e _ K2X / -Vl + sin^ + S2 sin0 2 \ {mx> Q 
y -s 2 v 1 + sin #2 - sin# 2 y 

Here, s 3 = sgn(E — Vj), hk y = s\y/\E - V±\"f/v 2 sin#i, 
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y/\E- V 3 \~f/v 2 Jl + sm 2 9 ] , j G {1,2} refers to the 
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FIG. 3: Dependence of the probabilities \a\ 2 and \b\ 2 of the 
evanescent waves on a at 9\ — n/4. 



region with Vj, and the propagation angle 9 2 of the 
transmitted wave is governed by the conservation of p y , 
siy/\E - Vi \ sin^i = s 2 yJ\E - V 2 \ sin# 2 . The coefficients 
n, t\, a, and b are determined by the continuity of 
i>K,k {x) and dipK,k [x)/dx at x = 0. We introduce a 
parameter a = (V 2 — V\)/(E — V\). the ratio of step 
height and the kinetic energy of the incident wave. The 
following effects depend only on a and 9\, regardless of 
other details such as the type (n-n or p-p) of the step. 

Phase jump and Berry phase. — We examine the de- 
pendence of the reflection amplitude n = |ri|e i6,r i on a 
and #1 G [0, 7r/2]; r\ is an even function of 0\. We first 
discuss the case of \a\ <C 1, and then that of finite a. 

For I a I <C 1, the reflection phase ri shows an abrupt 
jump of 7T at 61 — 7r/4; sec Fig. [5J The phase jump is 
accompanied by a reflection zero, |ri| 2 = 0. To see the 
behavior, we derive the expression of n near 9\ = 7r/4, 



cos26'i 

fi ~ ^ — a 

4cos 2 6»i 



for \a\ < 1 and 6 X ~ tt/4. (3) 



The 7T jump occurs, irrespective of the sign of a. And, 
for \a\ -C 1, the evanescent waves are ignorable (|a| 2 , 
|6| 2 ~ 0), as shown for 6\ = tt/4 in Fig. the same 
occurs for other values of 0\. Below, we attribute the tt 
jump to the Berry phase associated with the pseudospin. 

According to the Hamiltonian |T]), the pseudospin cou- 
ples with the vector q such that it is parallel (antiparallcl) 
to q in the n-n (p-p) step. We follow the change of q, to 
see the change of the pseudospin in the reflection. The 
vectors q of an incident wave with 6\ = ir/A — S (where 
(5 is a small positive angle) and the wave 4" / formed by the 
reflection of at the step are drawn in Fig. |H One can 
assign a clockwise path T to the change of q; the choice 
between clockwise and counterclockwise does not mat- 
ter. We also follow the change in the reversal-symmetric 
process Tq from 9^/ to 8*,. This process is also a solu- 
tion of Eq. ([1]) with the same energy, since the evanescent 
waves can be ignored and thus [Hk, 0] = 0. Notice that 
the incidence angle of 0$/ is 7r/4 + 6, and that F and 
Tq are reflection-symmetric about q y -axis. As &\ — > 7r/4 
(<5 — > 0), the spatial propagation in the process T becomes 
identical to that in Tq. However, at 0\ = tt/4, the change 
of the pseudospin in the T process differs from that in Tq 
by 2tt rotation, i.e., the difference T — Tq forms a loop 
encircling once the origin of the q space. The resulting 
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FIG. 4: (Color online) Schematic view of the change of <f in the 
reflection process for \a\ <C 1. This view also shows the change 
of pseudospin, since pseudospin is parallel (antiparallel) to q 
in the n-n (p-p) step. The change from an incident wave "Si 
(with 6\ < 7r/4) to its reflected wave ^ f is represented by 
r. The change of the reversal-symmetric process from O*]// 
(which incidence angle > 7r/4) to G^i is denoted by Tq. At 
9 1 - 7r/4, *i (*/) has the same q as 9*/ (6*;)> but the 
change of the pseudospin in the process F differs from that 
in Fq by 2n rotation, resulting in the abrupt tt jump (Berry 
phase tt in the q space) of the reflection phase 9 ri . 



Berry phase tt gives rise to the abrupt jump in 9 ri . 

We next discuss the case of finite a. As \a\ increases, 
the abrupt jump of 9 ri at 9\ = tt/A becomes gradual 
[Fig. [2J. This is due to the evanescent waves localized 
at the step boundary x = 0, which become to affect the 
reflection as |a| increases [Fig. [3], and break the rever- 
sal symmetry 0; for example, a wave ©| ± inj,k y ) oc 
\k y , -iziKj), reversal-symmetric to an evanescent wave 
| ±iKj 7 ky), is physically meaningless, as it diverges in 
y direction. For 9\ ~ 7r/4 and small |cv|, we derive 



dd r 



-(\/3a/12) 



- tt/A + a/4) 2 + (V3a/12) : 



+ 0(a 3 ), (4) 



which shows that the (gradual) jump of 9 ri occurs around 
a shifted angle 9\ = (tt — a) /4 within \/3|a|/12. For 
M ^ 1) @ri increases only gradually. Any bipolar step 
does not show the phase jump, as it has |a| > 1. 

Interferometry. — We propose an interferometry setup 
for the detection of the phase jump in 8 ri [Fig. [ljb)]. 
In addition to the step, it has a potential barrier V c , 
which shows transmission resonance only around some 
incidence angle 6 C [Fig. [5fa)]. This filtering or collima- 
tion Q is used in our setup, to see the dependence of 9 ri 
on the incidence angle (now on the collimation angle 9 C ). 

We first discuss the collimation. Figure [5ja) shows 
the dependence of the transmission probability T = \t\ 2 
through the setup on the incidence angle 9 unc of a plane 
wave. Here, t is calculated in the same way as for the 
step, and we choose (3 = (V c — V% )/{E — V\) > 2, for 
which there is no total reflection by the barrier for any 
#unc- The collimation angle 9 C , at which T shows a 
maximum value for a given (3, is governed by the res- 
onance condition, 2k Cx d c + ipo = lixn (n is an integer), 

where k c „ 
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FIG. 5: (Color online) (a) Transmission T = |t| 2 of a plane 
wave with incidence angle 9 unc through the setup in Fig. H^b) , 
for different /3's. T is an even function of 9 unc . (b) Conduc- 
tance G as a function of d s . In (a) and (b), E — Vi = 0.047, 
a = —0.5, and d c = 60io are chosen, while d s = 60io only in 
(a). Here, Go = e 2 W/(2n hd c ) and l = hv/j w 1.7 nm. 



Schn\/\/3 — 1| — sin 9 C is the i-axis wave 



vector inside the barrier, k- m = ^f\E — Vi\^/{hv) 2 is the 
wave vector of the incident wave, d c is the barrier width, 
s c = sgn(E — V c ), and tpo is the reflection phase at the 
barrier boundaries; in our parameter range, the depen- 
dence of (po on (3 and 9 unc is ignorablc. The collimation 
is mainly done in the barrier, and negligibly affected by 
the step, as |r*i | is small. One can tune 9 C by changing (3. 

Due to the collimation, a wave incoming into the step 
has incidence angle around 9 C . Its two propagation paths 
in x € [0, d s ], one with direct transmission and the other 
with reflection once at x = d s [Fig. |Tfb)] , result in the 
interference pattern of cos(2ki :c d s +9 ri +9 r2 ) as a function 
of d s , where hk-\ , = s\ki n cos9 c and 9 r2 is the reflection 
phase at x = [24|; for (3 > 2, 9 r2 changes only gradually. 
From the period A = 7r/fci x and the phase shift of the 
pattern, one can identify 9 C and the 7r jump of 9 ri . 

Based on the above idea, we analyze the conductance 
G = (4e 2 /h)(W/2ir) J dk y T(k y ) through the setup at 
zero temperature and zero bias, where the factor 4 re- 
flects the spin and valley degeneracy, the integral is done 
over all the incoming waves with the same energy E but 
different 9 unc , and W is the transverse width of the setup. 
For different /3's and thus c 's, the interference pattern in 
G is shown in Fig.[5jb), as a function of d s . The period 
of the pattern agrees very well with A = Tr/ki^, , therefore, 
from the period, one can identify the collimation angle 9 C 
in experiments, provided that fc; n is known. And, to see 
the tt jump of 9 ri , we redraw G in Fig. [SI as a function 
of d s /X; the jump is not clearly shown in Fig. E^b), since 
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FIG. 6: (Color online) The same as in Fig. Ob), but as a 
function of d a /\. 



A varies with (3. It shows the phase shift by tt around 
(3 = 4.05, where 9 C = I.IItt/A; the value l.ll7r/4 devi- 
ates from 7r/4, and the deviation roughly agrees with the 
deviation a/4 for small a [see Eq. (QJ]. In this way, by 
tuning P and thus 9 C , one can observe the tt jump of 9 ri . 

Discussion. — We include the band gap effect A(x) = 
Aj due to the external gates creating Vi. It is described 
by Hamiltonian ~ H K + A(x)a z /2 for A <C 7, 
and the argument in Fig. 0] has to be modified since 



^ for A ^ 0. According to Ref. (15J, A, 
can be expressed, using typical experimental parame- 
ters, as A ?; ~ -2&Vi, where &(e ?; ) = 0.5(1 + |e,|)/[l + 
l e i| + e i ~ 0.5 In |ei|] and ej = Vi/j. For small a [~ 
(ei — 62) /ei], we approximately use constant £j=i,2 = £, 
since £i(ei) varies slowly enough; we later consider the 
dependence of £j on ej. And, in the limit of zero bias 
and zero temperature, E ~ 0. In this regime of our 
interest, one can find the unitary pseudospin-rotation 
U satisfying UH^U f = (v 2 /j)a K ■ Q + V(x), where 
Q = (q/\q\)^q 2 +j 2 A , f/(2v 2 ) 2 = q/y/l^ for all i's. 
Notice that UH^W has the same form as Hk and that 
[J7H^C7t,e] = 0. Then the argument in Fig. H is ap- 
plicable for U^ij and QU^fij in the Q space, thus, the 
abrupt tt jump in 9 ri is maintained and can be detected as 
in Figs.[5]and[6l Eqs. ©-© and Figs. HO are not altered, 
and the dependence of 9 C on (3 is modified slightly. These 
features persist when we go beyond the approximation 
made above. For this case of £1 7^ £2, Eq. (|3|) is shifted 
as ri — > r\ + Sri, where Sri — — iasin20i/(4cos 2 9\) 
and a = Ai/[2{E - Vi)} - A 2 /[2{E - V 2 )}. From 
we find that Sri is ignorable (|a/a| < 0.1) in reasonable 
ranges of < 0.2, |otj < 0.5, and \E/V 2 \ < 0.2. This 
estimation of Sri will be modified only slightly when two 
(top and bottom) gates are used to create the step with 
small Vi —V2; in this general case, the above expression 
of £i(ei) may be altered, but the expression of Sri is still 
valid. Thus the tt jump in 6 ri is detectable in the pres- 
ence of the gap. We emphasize the central role of the Q 



space as the parameter space for the Berry phase tt. 

We compare our result with the Klein effect in mono- 
layer graphene. In the monolayer, it was predicted 11 1 
and observed [3| in an interferometry that a sign change 
(tt phase jump) occurs in the back-reflection amplitude. 
Contrary to our bilayer case, (i) its origin is the Berry 
phase tt in the p space, (ii) it occurs at zero incidence 
angle, thus an external magnetic field may be required 
to detect it, and (iii) it occurs in a bipolar junction. 

In summary, we find the abrupt jump tt of the reflection 
phase at a monopolar potential step in Bernal-stackcd 
bilayer graphene. The jump is the manifestation of the 
reversal symmetry and the Berry phase tt in the Q (or 
q) space, and robust against the band gap opening. We 
propose the setup for the detection of the jump, in which 
collimated waves are generated, tuned, and identified. 

This work is supported by NRF (2009-0078437). 
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